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Context. SGR 1806-20 is the soft gamma-ray repeater that exhibited the highest activity since its discovery, including a giant flare in 
2004, December 27. Previous studies of this source - which is likely a magnetar - showed that it was associated with a massive star 
cluster embedded into a gas and dust cloud. Moreover, several stars in the cluster are peculiar hypergiants - luminous blue variable 
(LBV) and Wolf-Rayet stars - exhibiting strong and likely dusty stellar winds. 

Aims. We aimed at studying the mid-infrared emission of the stars associated with the same cluster as SGR 1806-20, to detect 
variations that could be due to the high-energy activity of the magnetar through interaction with the dust. We also aimed at studying 
q , the morphology of the cloud close to the centre of the cluster. 

j_( . Methods. We carried out mid-infrared observations of SGR 1806-20 and its environment - with the highest spatial resolution in 

this domain to date - using ESO/VISIR in 2005 and 2006, and we retrieved 5p;?zer/IRAC-MIPS archival data of the same field. We 
performed broadband photometry of three stars - LBV 1806-20, a WC9 and an O/B supergiant - on our VISIR images, as well as 
on the IRAC data. We then built and fitted their broadband spectral energy distributions with a combination of two absorbed black 
bodies, representing their stellar components, as well as a possible mid-infrared excess, in order to derive their physical parameters. 
Results. We show that LBV 1806-20 and the WC9 star exhibit a mid-infrared excess, likely because of the presence of circumstellar 
dust related to their winds. We also show that only LBV 1806-20 had a variable flux over a period of two years, variability which 
is due to its LBV nature rather than to a heating of the gas and dust cloud by the high-energy emission of SGR 1806-20. Finally, 
differences in the intrinsic absorptions of the three stars show an inhomogeneous structure of the density of the gas and dust cloud in 
• the massive star cluster. 
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1. Introduction IWoods & Tho mpson 2006, for recent reviews). It is now com- 
monly accepted that SGRs and AXPs are magnetars, i.e. neutron 

Soft gamma-ray repeaters (SGRs) represent a small group - four stars whofje hi h _ e emission (persistent and transient) is 

£ i known objects, three in the Galactic centre and one in the Large powered by me decay of a verv strong dipole-like magnet ic field 

£ . Magellanic Cloud - of highly-magnetized exotic pulsars present- (Duncan & xhompson 1992; Thompson & Duncan 1995i 1996) . 

ing no evidence of binarity. They are characterised by a spin pe- , „„^„„. . „„„ . 

. 5 , , . ,. a it a ■* / ft v SGR 1806-20 is the SGR that ex hibited the highes t activ - 

nod clustered in the range 6-12 s, and a persistent soft X-ray ,. rr — : 71 1 . r " , ir;^o-^ 

• • , , n 35 ,^36 -1 to. • t 11 a ity since its discovery (Atteiaetal. 1987; Laros et al. 1987). 

emission around 10 - 10 erg s 1 . Their spectra are well-fitted T J . , , , J , ' , — : — 11 " ' ., . — , 

, .. c , , , , j n rr n c n /ci i?\ a It is located towards the Galactic centre, exhibits 7.47 s pul- 

by a combination of a black body (kT~0. 5-0.6 kev) and a power- . rr nr— tti. , „„,„ n 1 . 

, . ~, c , j , ..... , . j. . , . sations (Kouvehotou et al. 1998) and an 8.3 10 s s spin- 
law tail (r ~ 2.5), and they also exhibit short and intense bursts , * — : " . , . , , .. . . 

t . c c , , i .„• , , c c . au a down rate. It entered in an active phase in 2003 and exhibited 

(duration of a few hundred milliseconds) of soft y-rays and hard , , „. „ . „ , 

v /--.a u iu u a c \u c a. bursts up to late 2004. Finally, a giant flare (corresponding to 

X-rays. Giant flares have also been observed for three of them: , nA n f n A X _, „ ,. ■" b , „ , .K. , ^ 

c/^r> nc-i^ <lc cpd tont on a c/-tj mm 1 a r> *i 10 _ 10 erg s 1 for a distance between 8 and 15 kpc) oc- 

SGR 0526-66, SGR 1806-20 and SGR 1900+14. Recently, per- , . „„„„°. . ,,, . ^ , l u^od , j 

•• uu a i i a c/-^r. lonron curred in 2004 December 27 (Hurley et al. 2005), and a remnant 

sistent hard X-ray emission has been detected on SGR 1806-20 .. „ , „ , * , — . r , 

El . . . . . i l Unng k radio afterglow allowed an accurate localisation of the source 

( |Mereghettietal. ||2QQ5|K feaenS alJ l2005h . 

As for anomalous X-ray pulsars (AXPs) - proposed to be 

strongly related to SGRs due to their similar characteristics LBV 1806-20, a luminous blue variable (LBV), was long be- 

( Gavriil et al. 2002) - the likely source of their radiative emission p£ve& to be its near-infrared_(NIR) counterpart (Kul karni et al. | 

is not accretion like for other neutron st ars but their very strong but | Kaplan et al. | (|2002j) showed that the SGR was 12" 

magnetic field (10 14 - 10 15 G) (see e.g. lMereghettilt^l2rjof awa y from the LBV usln 8 Chandra observations. The accurate 

localisation of the source, as well as an exhaustive monitoring fi- 

Send offprint requests to: F. Rahoui nally led to t he detection of the NIR counterpart of SGR 1 806-20 

* Based on observations carried out at the European Southern (Israel et al.||2005t |Kosugi et al.||2005|) by observing a variabil- 

Observatory under programmes ID 075.D-0773 and 077.D-0721 (P.I. ity correlated to the high-energy emission, as for AXPs, which 

S. Chaty). strengthens the connection between both objects. 
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The environment of SGR 1806-20 is very young and dusty, 
and it is then relevant to observe it at mid-infrared (MIR) wave- 
lengths and search f or variabiliti e s corre lated t o the high-energy 
emissi on. Indeed, iFuchs et ail d 19991) and lEikenberrv et al.l 
(1200 11) showed that it was associated with a massive star cluster, 
along with LBV 1806-20. Moreover, this cluster is located 
inside the radio nebu la G10.0-0.3, powere d by the very strong 
winds of the LBV dGaensler et alj 12001 1) . Finally, this radio 
nebul a itself is associated with the giant molecular complex 
W31 (ICorbel et al]|1997h . ICorbel & Eikenberrvl d2004l) showed 
that W3 1 was resolved into two components, one at d ~ 4 kpc 
withA v ~ 15, and another one atd~ 15.1 kpc with A v = 37 + 3, 
and that LBV 1806-20 (and consequently the massive star 
cluster, and SGR 1806-20) belonged to the further component 
of W3, sug g esting a distance of 15. 1*H kpc. Nevertheless, 
Bibb v et ail (|2008) recently performed high-resolution near- 
infrared spectroscopy of O/B and Wolf-Rayet stars in the cluster 
to derive their accurate spectral classification. This allowed 
them - using both synthetic photometry and isochrone fitting - 
to derive a distance modulus of 14.7 ± 0.35 mag for the cluster, 
which corresponds to a distance of 8.7+J ■? kpc. 

In this paper, we report MIR observations - with the highest 
spatial resolution in this domain to date - of SGR 1806-20 and its 
environment, carried out at ESO/VLT with VISIR in June 2005 
and 2006. We also report archival data obtained at three different 
epochs between 2004 and 2006 with IRAC and MIPS mounted 
on Spitzer. Our goal was to detect any MIR variability which 
could be due to the high-energy activity of the SGR, as well as 
to study the morphology of the gas and dust cloud in which all 
sources are embedded. 

In section 2, we describe these observations and their analy- 
sis. In section 3, we present the broadband spectral energy dis- 
tributions (SEDs) we built for three sources using the VISIR and 
IRAC fluxes obtained in 2005 and 2006. For each year, there 
were no more than three months between the VISIR and the 
IRAC observations, which allowed us to fit these SEDs to de- 
rive the physical parameters of the stars with contemporaneous 
data from 3.6 /an to 11.25 /mi. In section 4, we analyse the re- 
sults and discuss the MIR variability, the intrinsic absorption, 
and the distance of these stars, as well as the MIR emission of 
SGR 1806-20. 



2. Observations and data reduction 

2.1. ESO/VISIR data 

The MIR observations were carried out on 2005 June 2 0-22 
and 2006 June 29-30 using VISIR dLagage et all 12001 . the 
ESO/VLT mid-infrared imager and spectrograph, composed of 
an imager and a long-slit spectrometer covering several filters in 
N and Q bands and mounted on Unit 3 of the VLT (Melipal). 
The standard "chopping and nodding" MIR observational tech- 
nique was used to suppress the background dominating at these 
wavelengths. Secondary mirror-chopping was performed in the 
north-south direction with an amplitude of 16" at a frequency of 
0.25 Hz. Nodding technique, needed to compensate for chopping 
residuals, was chosen as parallel to the chopping and applied us- 
ing telescope offsets of 16". Because of the high thermal MIR 
background for ground-based observations, the detector integra- 
tion time was set to 16 ms. 

We performed broadband photometry in three filters, 
PAH1 (.1=8.59+0.42 /mi), PAH2 (^=11.25+0.59 /an), and 
Q2 (A= 18.72+0.88 /mi), using the large field in all bands 



(32'.'5x32'.'5 and 07127 plate scale per pixel) in June 2005. In 
June 2006, we used the large field in PAH1 and PAH2 and the 
small one (1972x1972 and 07075 plate scale per pixel) in Q2. 
All the observations were bracketed with standard star observa- 
tions for flux calibration and PSF determination. The weather 
conditions were good and stable during the observations. 

Raw data were reduced using the IDL reduction package 
written by E. Pantin. The elementary images were co-added in 
real-time to obtain chopping-corrected data, then the different 
nodding positions were combined to form the final image. The 
VISIR detector is affected by stripes randomly triggered by some 
abnormal high-gain pixels. A dedicated destriping method was 
developed (Pantin 2008, in prep.) to suppress them. 

The filtered reduced images of 2006 are displayed in 
Fig 1. These images were cleaned using the multiresolu- 
tio n software package writ ten by JL. Starclfl and presented 
in IStarck & Murtaghl {1998), which performs background and 
noise modelling, as well as noise subtraction using multires- 
olution tools. Nevertheless, we performed photometry on the 
unfiltered images using aperture photometry, and the fluxes 
in all bands are lis t ed in Table 1. We kept the s ame labels 
as lEikenberrv et all (120041 EML04 hereafter) and iFiger et ail 
(2005]) for all the detected stars. A is therefore LBV 1806-20, 
B is a WC9 star, C is an O/B supergiant, and 5 is a Red Giant. 
Finally, G is a hot and dense part of the gas and dust cloud which 
is the likely birth site of the massive star cluster, and is only vis- 
ible in the MIR domain. On June 2005, the exposure time was 
600 s in all bands while it was 1200 s in PAH1 and PAH2 and 
2400 s in Q2 in June 2006, which explains why we did not de- 
tect neither C in PAH1 nor G in PAH2 in June 2005, their fluxes 
being too low compared to the sensitivity we were able to reach 
(see Table 1). A three-color image (made with PAH1, PAH2, and 
Q2 images) is also displayed in Fig 2 and we clearly see that the 
gas and dust cloud emits mostly beyond 1 1 /mi. 

2.2. Spitzer/IRAC-MIPS archival data 

We searched for MIR public data of the star cluster in the 
archives of the Spitzer space telescope. We found photomet- 
ric data taken at three different epochs, 2004 October 29, 2005 
September 27, and 2006 April 10, in the Spitzer's, Galactic 
Legacy Infrared Mid -Plane Survey Extraordinaire (GLIMPSE, 
Benjamin et al. 2003), survey of the Galactic plane (\b\ < 1° 
and between /=10° and /=65° on both sides of the Galactic 
centre), using the IRAC camera in four bands, 3.6+0.745 /mi, 
4.5±1.023 /mi, 5.8+1.450 /an and 8.0+2.857 /an. IRAC plate 
scale is 172 per pixel and the field of view is 5'2x5f2. We also 
found archival MIPS data of the same field at 24 /an. 

We performed photometry of the stars A, B and C we had 
detected with VISIR on the post-Basic Calibration Data (post- 
BCD) using the software StarFinder, part of the Scisoft package 
from ESO, well-suited for point-source photometry in crowded 
fields. Post-BCD data are raw data on which the Spitzer pipeline 
performs dark subtraction, multiplexer bleed correction, detec- 
tor linearization, flat-fielding, cosmic ray detection, flux calibra- 
tion, pointing refinement, mosaicking, and coaddition. The im- 
ages taken in September 2005 are displayed in Fig 3, and fluxes 
for each epoch are listed in Table 2. Moreover, Fig 4 displays the 
MIPS image of the same field at 24 /an. 



1 http://thames.cs.rhul.ac.uk/~multires/mr4-software/index.html 
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Fig. 1. PAH1 (left, a.' 127 plate scale per pixel), PAH2 (centre, a.' 127 plate scale per pixel), and Q2 (right, a/075 plate scale per 
pixel) VISIR images of the environment of SGR 1 806-20 (white X point) obtained in J une 2006. A is LBV 1806-20, B is a WC9 
star, C is an O/B supergiant, and 5 is a Red Giant (Eikenberry et al. 2004]; iFiger et a"Dl2005l) . An ellipse pattern surrounds G, the 
hottest part of the gas and dust cloud in which all stars are embedded. We detected it in PAH2 and in Q2. 
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Fig. 2. Three-color image of the environment of SGR 1806-20 (white X point) done with images of June 2006; PAH1 is in blue, 
PAH2 in green, and Q2 in red. Q2 image was rescaled to a Of.' 121 plate scale. The cloud detected with VISIR appears to be splitted 
into several components, which shows that its temperature is not homogeneous, the cloud being hotter in the centre of the massive 
star cluster. 
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Fig. 3. IRAC 3.6 /urn (top left), 4.5 fim (top right), 5.8 /vm (bottom left) and 8.0 //m (bottom right) images of the environment 
of SGR 1806-20 (white X point). The field of view is about 79"x71" and the plate scale is 1'.'2 per pixel. We point out that the 
position of SGR 1806-20 is 2'.'95 away from the point source located i south-east of the white X point, which - considering the (X'5 
astrometric accuracy of the Spitzer telescope - excludes an association with SGR 1806-20. 
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Fig. 4. MIPS 24 /im image of the environment of SGR 1806-20 (black X point). We also superimposed the MIPS contour plot to 
stress the differences in luminosity. We clearly see the gas and dust cloud in which all the stars are embedded, and which is likely 
the birth site of the massive star cluster. It appears that LBV 1806-20 (A) is located in a less dense and colder zone of the cloud than 
the WC9 star (B) and the O/B supergiant star (C), which explains that its intrinsic absorption is lower. Moreover, the central part of 
the cloud - which appears saturated in this image as it is too hot with respect to the MIPS level of saturation - corresponds to the 
region detected with VISIR. The field of view is 94"x67" and the plate scale is 2'.'5 per pixel. 

Table 1. Summary of VISIR observations of the sources in the environment of SGR 1806-20. We give their name and their MIR 
fluxes (in Jy) in the PAH1 (8.59 //m), PAH2 (11.25 /mi), and Q2 (18.72 fim) filters. When we did not detect the source, we give 
its upper limit. For the dust cloud, we give the upper limit for an extended source covering 450 pixels in PAH1 and PAH2, as this 
corresponds to the area covered by the dust cloud when we detected it in PAH2 in June 2006. Uncertainties are given at lcr. 



Sources PAH1 


PAH2 


Q2 


June 2005 


LBV 1806-20 (A) 0.133 ± 0.004 


0.089 ± 0.005 


< 0.058 


WC9 (B) 0.077 ± 0.003 


0.057 ± 0.004 


< 0.058 


Supergiant (C) < 0.0126 


< 0.0149 


< 0.058 


Red Giant (5) 0.020 ± 0.002 


< 0.0149 


< 0.058 


Dust (G) < 0.091 


< 0.111 


2.60 ± 0.20 


June 2006 


LBV 1806-20 (A) 0.103 ± 0.004 


0.067 ± 0.004 


< 0.051 


WC9 (B) 0.076 ± 0.003 


0.053 ± 0.003 


< 0.051 


Supergiant (C) 0.0 1 2 ± 0.002 


< 0.0124 


< 0.051 


Red Giant (5) 0.019 + 0.002 


< 0.0124 


< 0.051 


Dust(G) < 0.112 


0.078 ± 0.007 


2.82 + 0.14 



3. Absorption and SEDs 

3.1. Absorption 

The absorption at wavelength A, A,\, is a crucial parameter to fit 
the SEDs, especially in the MIR. Indeed, non accurate values 
can lead to a false assessment of the dust emission. The visible 
absorption A v was a free parameter of the fits. An accurate in- 
terstellar absorption law - i.e. the wavelength dependence of the 
j 1 ratio in the line of sight - was then needed to properly fit the 
SEDs. 

From 1.25 /um to 8.0 / / m, we used the analytical expression 
given in Indeb etouw et al.l (120051) . They derived it from the mea- 
surements of the mean values of the color excess ratios /~. J 
from the color distributions of observed stars in the direction of 
the Galactic centre. They used archival data from the GLIMPSE 
catalogue, which is relevant in our case. Beyond 8.0 /mi, where 
absorption is dominated by the silicate feat ures at 9.7 um an d 
18.0 /mi, we used the extinction law from iLutz et al.l ([1996), 



which includes the interstellar silicate absorption at 9.7 /mi. The 
values we used in each band are listed in Table 3. 



3.2. SEDs 

With all the observational and archival data from NIR to MIR 
wavelengths, we built the 2005 and 2006 SEDs of these sources. 
We fitted them (x 1 minimisation) using a combination of two ab- 
sorbed black bodies, one representing the companion star emis- 
sion and one representing a MIR excess, if there was any: 



AF(A) 



2nhcl lQ -0AA A 
D t 2 A 4 



R* 



Ru 2 



he 
,AkT, 



he 
,AkT D 



in WirT 2 (l) 



1 e AkT o - 1- 

We added to the data uncertainties systematic errors as follow- 
ing: 

- a 2% systematic error in each IRAC band as given in the 
IRAC manuafl 



http://ssc.spitzer.caltech.edU/documents/som/som8.0.irac.pdf 
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Table 2. Summary of IRAC observations of the sources in the environment of SGR 1806-20. We give their name and their MIR 
fluxes (in Jy). Uncertainties are given at lcr. 



Sources 


3.6/um 


4.; 


i yum 


5.8/im 


&.0imi 


October 2004 


LBV 1806-20 (A) 


0.306 


± 


0.010 


0.316 


± 0.008 


0.283 ± 


0.008 


0.175 


± 


0.005 


WC9 (B) 


0.165 


± 


0.005 


0.208 


± 0.006 


0.218 ± 


0.007 


0.141 


± 


0.003 


Supergiant (C) 


0.041 


± 


0.002 


0.036 


± 0.002 


0.028 ± 


0.002 


0.027 


± 


0.002 


September 2005 


LBV 1806-20 (A) 


0.337 


± 


0.015 


0.332 


±0.016 


0.299 ± 


0.014 


0.199 


± 


0.007 


WC9 (B) 


0.172 


± 


0.009 


0.218 


±0.011 


0.224 ± 


0.011 


0.140 


± 


0.005 


Supergiant (C) 


0.037 


+ 


0.003 


0.037 


± 0.003 


0.031 ± 


0.003 


0.027 


± 


0.001 












April 


2006 










LBV 1806-20 (A) 


0.297 


± 


0.010 


0.304 


±0.001 


0.272 ± 


0.011 


0.174 


± 


0.005 


WC9 (B) 


0.170 


± 


0.006 


0.207 


± 0.007 


0.225 ± 


0.009 


0.140 


± 


0.006 


Supergiant (C) 


0.041 


± 


0.003 


0.038 


± 0.002 


0.030 ± 


0.003 


0.025 


± 


0.002 



- comparing the variations of the flux calibration values ob- 
tained from the standard stars with VISIR during our ob- 
servation nights, we figured out that systematic errors were 
about 5% at 10 /im and 10% at 20 fim. 

The free parameters of the fits were the absorption in the V band 
A v , the companion star black body temperature T„ and its radius- 
to-distance ratio as well as the MIR excess spherical compo- 
nent black body temperature and radius To and Rr,. The best- 
fitting parameters for individual sources, as well as the corre- 
sponding^ 2 are listed in Table 4, and fitted SEDs of objects A, 
B and C are displayed in Fig 5. 

4. Results and discussion 

4.1. The variability of LBV 1806-20 

The fluxes listed in Tables 1 and 2 show that LBV 1806-20 is 
variable in the MIR, in all bands from 3.6 //m to 11.25 fim. 
Indeed, The IRAC fluxes substantially increased from October 

2004 to September 2005, and then decreased from September 

2005 to April 2006 to reach a similar level as in October 2004. 
Concerning the VISIR fluxes, we unfortunately do not have any 
data taken in 2004. Nevertheless, the same decrease in the MIR 
flux is observed from June 2005 to June 2006 in both PAH1 and 
PAH2. 

In a previous spectroscopic and photometric study of 
LBV 1806-20, EML04 confirmed the likely LBV nature of this 
star (van Kerkwijk et al., 1995) with spectral classes 09-B2, and 
showed that there were strong anticorrelated variations of the 
equivalent widths of He i 2.1 12 //m and Bry lines, which led to 
correlated variations of the star temperature and the number of 
ionisi ng photons. Moreov er, comparison with previous photom- 
etry dKulkarni et aljfl995l) showed that the variation of temper- 
ature did not result in the variation of the K band magnitude of 
the star. As a possible explanation, they proposed either an an- 
ticorrelated variation of the star temperature and radius, which 
is typical of some LBV stars dHumphrevs & Davidson! IT994t 
Morr is et"al][T9 96), or an increase in the stellar wind density, 
resulting in more absorbed ultraviolet continuum photons and 
consequently, more emission lines. 

To investigate the origin of the MIR flux variation of 
LBV 1806-20, and to find out whether it was related to one of 
the previously described mechanisms or to the heating of the gas 
and dust cloud by the high-energy emission of SGR 1806-20, we 
built its broadband SED using 2005 and 2006 IRAC and VISIR 



data - as well as the NIR magnitudes from EML04, the NIR flux 
being marginally variable - and fitted it. The best-fitting param- 
eters are listed in Table 4. 

The presence o f warm circumstellar dust around LBV stars is 
common (see e.g. McGregor et al. 1988, about AG Car), as they 
are often associated with ejected dusty nebulae, and a compo- 
nent of warm dust is necessary to explain the huge MIR excess 
LBV 1806-20 exhibits. 

Moreover, the fits suggest a change from a Bl (T* = 
21500 K) to a B3 (r» = 17600 K) spectral type, and at least 
a 20% increase in the stellar radius (R t = 3.08 x 10~ 10 x D t . 
to R* = 3.68 x 10~ 10 x £>„) from 2005 to 2006. Meanwhile, 
the dust temperature decreased by about 10% (7b = 1130 K to 
To = 1020 K) while the dust radius (decrease of about 7.5% 
fromtfo = 25.87 x 10 -10 x D t to R D = 23.92 x 10~ 10 x D t ) 
was roughly constant, considering the uncertainties on the star 
radius and distance (about 17%). 

Under the assumption that a dust grain is a perfect black 
body - i.e. it fully absorbs the received flux - and that the star 
emission is not absorbed before reaching the dust, the equation 
of thermodynamical equilibrium between a dust grain and the 
star can be written: 



Rd 



4x Q(v)B(T D ,v)dv 



(2) 



where B(Tu, v) is the dust grain black body emission at the 
frequency v, and Q(v) the dust grain emissivity, approximated 
as Q(v) - Qti^, where Qq is a constant and 1 < n < 2 
dDraine & Leefl984l) . 

The expected grain temperature can therefore be written: 



1 



-(-) 

6060 A*/ r(4+n)£(4+n)\RD) 



-\ Tt 



(3) 



where h and k are the Planck and the Boltzmann constants, re- 
spectively, F the gamma function, and ( the Riemann zeta func- 
tion. 

By studying t he properties of the war m dust around the 
LBV HD168625. lRobberto & Herbstl dl998l) give On = 1.52 x 
10~ 8 x a - where a is the dust grain radius - and find a ~ 1 /mi, 
consistent with LBVs h aving large dust grains (a > 0.1 /mi) in 
thei r winds, as show n in Mitchell & Robinson ( 1986) for rj Car, 
and lMcGregoretal.1 d 1988b and lShore et all d 19961) for AG Car, 
respectively. 
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Table 3. Adopted ^ values. 



Filters J H Ks 3.6yum 


4.5 fim 5.8//m 


8.0yum 


8.59yum 


1 1.25 yum 


& 0.289 0.174 0.115 0.064 


0.054 0.047 


0.044 


0.060 


0.061 



Table 4. Summary of parameters we used to fit the SEDs of the sources. We give their name and then the parameters themselves: 
the extinction in the optical A v , the temperature T„ and the ratio of the star, and the temperature Tu and radius Ru (in R„ unit) of 
the dust component. We also add the reduced^- 2 we reach for each fit, as well as the uncertainties on the parameters. 



Sources 


A v 


r,[AT»](K) 


R, 
D, 




r D (K) 


R D (R.) 


^ 2 /dof 


LBV 1806-20 (A) (2005) 


28.6+^ 


21500 [17900 - 26200] 


3.08!° : *> x 10- 


•10 


i no +13 ° 


8 4 +1 


3.1/4 


LBV 1806-20 (A) (2006) 


29.0+1 1 


17600 [15000 - 21200] 


3.681°^ x 10" 


-id 


io2o:;s 


6 5 +07 


3.9/4 


WC9 (B) 


35.2+j| 


53000 [25000 - 71000] 


4.36^ x 10" 


-li 


1850% 


46 +8 ° 


6.8/4 


Supergiant (C) 


31.8!™ 


25700 [11600 - 40500] 


1.38+^ x io- 


-in 






6.0/5 



With Qq given in iRobberto & Herbstl d!998l) . a large dust 
grain size a ~ 0.5 fim, and n = 1.2 as given in 
iRouleau & Martini (1 19911) . we find T D ~ 1 143 K in 2005, and 
Tu ~ 1081 K in 2006, i.e. a decrease of about 6%, in good 
agreement with what we obtain in our best fits. This is therefore 
consistent with a MIR flux variation of LBV 1806-20 due to the 
anticorrelated variations of the star temperature and radius rather 
than to a heating or a cooling of the gas and dust cloud due to 
the high-energy activity of SGR 1806-20. 

4.2. Stars B and C 

B was classified as a WC9 Wolf-Rayet (EML04), and the fluxes 
listed in Table 1 and 2 show that its MIR emission was constant 
in all bands from 3.6 //m to 11.25 fim. We then fitted its SED 
using the MIR fluxes of 2006 and the NIR magnitudes found in 
EML04, and the best-fitting parameters are listed in Table 4. Our 
results suggest that B is extremely reddened and that it exhibits 
a very large MIR excess. The presence of such excess is typical 
of WC9 stars, whose C rich stellar winds allow the creation 
of circumstellar dust. Moreover, Fig 5 shows that beyond 
1.6 fim (H band), the star barely contribu tes to the emission, 
which is also common fo r WC9 stars (Ivan der Hucht et all 
Il996t ICrowther et al.ll2006l) . All the fits with a more balanced 
contribution of the stellar and the dust components (reduced 
X 1 > 3), or with a stellar component only (reduced x 2 ^ 100) 
failed to reproduce the SED. Moreover, we point out that 
although To = 1850K is high, it is still below the sublimation 
temperature of amorphou s carbon found around WC9 stars 
(-200 K, see for instance ILaor & Draindfl993l: iPreibisch et alj 
119931) . 

C was found to be either an O/B supergi ant or an hyper- 
giant through spectroscopy (EML04), then Fi ger et aTl ([2005 ) 
performed higher resolution spectroscopy and derived an O/BI 
spectral class, with n arrower lines than previously found. Their 
result is confirmed bv lBibbv et al.l d2008l) who derive a B1-B3I 
spectral type through high-resolution NIR spectroscopy. We un- 
fortunately did not detect C during our first VISIR run in June 
2005 as our exposure time was not sufficient, preventing us to 
perform a comparison with the flux of June 2006. Nevertheless, 
the IRAC fluxes listed in Table 2 point out towards a constant 
MIR flux from 2004 to 2006 in all bands. Using the NIR mag- 
nitudes given in EML04, as well as the IRAC and the VISIR 



fluxes of 2006, we fitted its SED and the best-fitting parameters 
are listed in Table 4. 

The best-fitting temperature is typical of an early-B su- 
pergiant star. Moreover, if we assume that the star C is asso- 
ciated with the cluster, it is then possible to get its radius by 
multiplying the best-fitting radius-to-distance ratio derived from 
our fit by the distan ce of the cluster. If we use a distance of 
15.1+11 kpc given in Corbel & Eikenberry (2004), we derive a 
radius R t = 92.4 + ^- 2 r Qi while we find R* = 53.2+^ R for 
a distance of 8.7+} * kpc dBibbv et al.ll2008l) . The latter is more 
consistent with the expected radius of a normal B 1-B3 super- 
giant s tar, therefore favouring the distance given in Bibb v et all 
d2008h . 

Finally, we would like to point out that the three-color image 
displayed in Fig 2 interestingly suggests that C is associated with 
the northern part of the gas and dust cloud, since the cloud emis- 
sion reappears at the supergiant position. If this were confirmed, 
it would mean that the cloud is heated by C. 

4.3. SGR 1806-20 

High-energy emission of magnetars is believed to be powered 
by the magnetic energy of their very strong magnetic field. 
Moreover, five AXPs and SGR 1806-20 have been detected in 
the optical bands and/or the NIR, and all of them were found to 
exhibit a variable NIR emission correlated to the X-ray emission. 
This optical/NIR emission is explained either as the non-thermal 
radiation by particles in the magnetosphere, or as the irradiation 
by the X-ray emission of a fossile passive dust disk around the 
neutron star. 

The second explanation recently received some credit as the 
AXP 4U 0142+61 h as been detected i n the MIR using IRAC 
at 4.5 fim and 8 fim (Wang et al. 2006). It is shown that its IR 
emission can be understood by the presence of an X-ray heated 
disk. This disk would have formed from the fallback material of 
a supernova, and would emit mostly in the IR consequently to 
the irradiation by the high-energy emission of the neutron star. 
Moreover, the authors derive an unabsorbed MIR to X-ray flux 
ratio V4 5 / , m f4 5 / , m /Fx ~ 3.4 x 10~ 4 , where Fx is the unabsorbed 
2-10keV X-ray flux. 

SGR 1806-20 was never detected in the MIR domain, nei- 
ther with IRAC nor with VISIR, as shown in the images dis- 
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Fig. 5. NIR to MIR fitted SEDs of LBV 1806-20, the WC9 star, 
and the O/B supergiant. For LBV 1806-20 and the WC9 star, 
we also add parts of their emission due to the stellar and the dust 
components. We used the MIR data from the VISIR observations 
we performed, the GLIMPSE survey, as well as the NIR data 
from EML04. The silicate absorption feature at 9.7 fan comes 
from the use of the Lutz's interstellar extinction law, in which it 
is included. 



played in this paper. Indeed, considering that the MIR emis- 
sion of SGR 1806-20 would be due to an irradiated fallback 
disk, and that the u nabsorbed MIR to X-ray flux ratio derived in 
Wan g et ail (|2006) is typical of dust disks around young neutron 
stars, we can expect the SGR 1806-20 absorbed flux at 4.5 jjm 



to be below 10 ply, one order of magnitude below the sensitivity 
of about 200 (dy at 5cr of the GLIMPSE survey. 

4.4. Absorption 

All the intrinsic absorptions derived from our fits show that stars 
B and C are likely more absorbed than LBV 1806-20, as they are 
embedded in a hotter and denser part of the associated cloud of 
dust and gas visible on our VISIR images in PAH2 and Q2 (see 
Fig 2). Indeed, Fig 4 displays a MIPS image of the environment 
of SGR 1806-20 at 24 pm with contour plots. The Spitzer sen- 
sitivity being far better, a larger part of the cloud is visible. We 
see that LBV 1806-20 (A) is located in a colder and less dense 
zone than the WC9 star (B) and the B supergiant star (C), which 
both are closer to the core of the cloud. Note that the image is 
saturated, and the saturated pixels correspond to the zone of the 
cloud detected with VISIR (taking into account the difference of 
plate scale per pixel in both images, 07075 for VISIR and 2'.'5 
for MIPS). B is then likely very close to the hottest and densest 
zone of the dust cloud, which is why it exhibits the highest in- 
trinsic absorption. We nevertheless point out that a limitation to 
this conclusion is that we only see the cloud integrated along the 
line of sight. 

5. Conclusions 

We reported the mid-infrared photometry of three stars and the 
gas and dust cloud associated with the same massive star clus- 
ter as SGR 1806-20, using ESO/VISIR and Spitze r/lRAC-MlPS 
data obtained at different epochs between 2004 and 2006. 

We show that LBV 1806-20 is the only object to exhibit a 
likely MIR variability, and that this flux variation is probably 
the consequence of its LBV nature rather than a heating of its 
circumstellar dust by the giant flare exhibited by SGR 1806-20. 
We also show that the stars in the central zone of the massive star 
cluster appear more absorbed, as they are closer to the hottest 
and densest part of the gas and dust cloud in which all of them 
are embedded. 

Finally, we recommend further high-sensitivity and long ex- 
posure MIR observations of SGR 1806-20 in order to try to de- 
tect it, and perhaps constrain the origin of the AXPs and SGRs 
optical/IR emission. 
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